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ABSTRACT 


The rapid growth of technology makes it possible to implement in immediate 
diagnosis for patients using image processing. By using morphological 
processing and adaptive thresholding method for segmentation of optic disc 
and optic cup, various sizes of retinal fundus images captured through fundus 
camera from online databases can be processed. This paper explains the use 
of color channel separation method for pre-processing to remove noise for 
better optic disc and optic cup segmentation. Noise removal will improve 


Keywords: image quality and in return help to increase segmentation standard. Then, 
Adaptive Thresholding morphological processing and adaptive thresholding method is used to 
| extract out optic disc and optic cup from fundus image. The proposed method 
Glaucoma is tested on two publicly available online databases: RIM-ONE and 
Morphological DRIONS-DB. On RIM-ONE database, the average PSNR value acquired is 
Segmentation 0.01891 and MSE is 65.62625. Meanwhile, for DRIONS-DB database, the 
best PSNR is 64.0928 and the MSE is 0.02647. In conclusion, the proposed 
method can successfully filter out any unwanted noise in the image and are 
able to help clearer optic disc and optic cup segmentation to be performed. 
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1, INTRODUCTION 

Glaucoma is one of the factors for irreversible blindness in the world. Visual field loss, optic disc 
cupping, and in most instances an increased intraocular pressure are the main features of glaucoma [1]. It is 
originated from increasing pressure called intraocular pressure (IOP). The pressure happens because of the 
leftover fluid that is gradually accumulated in the front part of the eye. The liquid cause increasing pressure 
within the eyeball and eventually lead to optic nerve damage. 

If immediate treatment was not delivered, it could lead to permanent eye damage and untreatable 
vision loss. Figure | shows glaucoma progression from normal vision until extreme glaucoma. 
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Figure 1. Normal vision progression into extreme glaucoma [2] 


However, not all high eye pressure caused glaucoma. Sometimes glaucoma can also happen to 
normal eye pressure as well. This depends on the level of pressure employed on the particular optic nerve. 
This pressure was usually measured using Goldmann Applanation Tonometry. Goldmann Applanation 
Tonometry was appraised as the main standard for IOP measurement [3]. 

Generally, there are two common glaucoma types which are primary open angle glaucoma (POAG) 
and angle closure glaucoma (ACG). About two thirds of glaucoma cases are POAG while the rest are ACG 
cases [1, 4-7]. Commonly, some well-developed health facility will provide free screening to detect 
glaucoma in patients [8]. This screening will be conducted using a fundus camera and performed by 
approved professionals and medical practices. 

For glaucoma infected eyes, it can be seen from the size of the optic cup and optic disc as illustrated 
in Figure 2. Optic disc is the path where ganglion cell axons exit and leave the eye. Optic cup is situated in 
front of the optic disc. It is shaped like a small crater. Conventionally, the size of optic cup is smaller than 
optic disc. Only in glaucoma infected eye, the optic cup size increases and grow larger. This is due to 
elevation of IOP. Besides, blood flow to the optic nerve is loss and nerve fibers slowly dying.The normal size 
for optic cup is 0.3 to 0.4 while glaucomatous eye have cup size of 0.7 to 0.9 [9]. 





Figure 2. Normal cup size versus glaucomatous cup size [10] 


The algorithm for this software that will be developed can perform the image processing techniques 
needed for the detection of glaucoma. There are two important steps completed. The steps are pre-processing 
phase and segmentation phase shown in Figure 3. 


Online Database 
RIM-ONE DRIONS-DB 


Pre-processing Phase 


Color Channel Separation 
Method 


Segmentation Phase 


Adaptive Thresholding and 
Morphology Method 





Figure 3. Block diagram of proposed system 
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2. RESEARCH METHOD 

The images used in this project were obtained from two online databases which are RIM-ONE [11] 
and DRIONS-DB [12]. RIM-ONE database consists of 169 Optic Nerve Head (ONH) images retrieved from 
169 full fundus images of different subjects. Figure 4 displayed examples of images from RIM-ONE. 
Domain experts have categorized the images into five different subsets: 
11 images of Ocular hypertension (OHT). 
14 images of Deep glaucoma. 
118 images of Normal eye (non-glaucomatous). 


12 images of Early glaucoma. 


14 images of Moderate glaucoma. 
(ii) (111) (iv) (v) 


(i) 


Figure 4. The examples of images from RIM-ONE database for (1) Ocular hypertension (OHT) (11) Deep 
glaucoma (111) Normal eye (non-glaucomatous) (iv) Early glaucoma and (v) Moderate glaucoma 


eno gp 





Meanwhile th next database which is DRIONS-DB demonstrates the segmentation of the optic 
nerve head in retinal fundus images. It was gathered by 2 medical experts. It has 2 ground truths for each set 
of 110 colour digital retinal images. Figure 5 illustrate the examples of images from this database. 


Figure 5. Samples of images from DRIONS-DB database 





2.1. Pre-Processing 

Pre-processing image is a common early method to obtain improved and enhanced image and to 
draw out useful information from it. Pre-processing usually involve images with the lowest level of 
abstraction. The aim of pre-processing is to remove noise, correct some of the degradation in image, 
suppresses unwanted distortions and amplifies image features important for further processing. 

For pre-processing, the method used is color channel separation method [13]. The fundus images 
obtained from databases are split into three color channels. Normally used color channels are red, green, and 
blue color channels (RGB). The RGB components usually regarded as three sub-sources of single color 
source [14]. 24-bit RGB image can give storage pixels with conventional brightness intensities of O until 255 
[15]. It first started with black color and the light emission increase as correspond to added colors. Figure 6 
shows examples of images obtained as the outcome of color channel separation method. 

One of the color channels is chosen and then processed further into grayscale. For this system, red 
color channel selected and undergoes the next process. 
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Retina Image Red Image 
Green Image Blue Image 





Figure 6. Fundus image after separated into three color channels 


2.2. Segmentation Phase 

The next process are segmenting and extracting important features from the grayscale image. 
obtained from pre-processing. Segmentation phase will separate an image into pixels with similar attributes. 
The pixels then accumulated into its distinguishable regions. This distinguishable regions need to have large 
relation to characteristic of interest for an accurate image interpretation. Typically for glaucoma, the 
important features that need to be extracted are blood vessels and optic disc. 

In this phase, the obtained grayscale image is converted into binary disc image. Binary disc image 
demonstrate raw disk images after conversion from grayscale image. 

After that, a further processed detailed image almost similar to binary disc image acquired using 
Morphological Processing method [16]. This image acquired 1s more refined that previous binary disc image 
existed as it improves and enhance image so that further processing can be made with ease. Generally, there 
are two formulas used in Morphological Processing method [16]. Each one represents dilation operation and 
erosion operation. Dilation operation is a type of operation that increase object in a binary image whether in 
terms of size or thickness meanwhile erosion operation is the operation that does the opposite of dilation 
method for example the operation reduce object’s size or thickness in binary image. These two operations are 
very important in increasing the accuracy of Morphological Processing method. It is expressed as: 


XB = X+b = {x+b:xEX&bEB} (1) 
Where X is any grayscale shape and B is symmetric structuring element. 

XOB = X-b = {z:(B+z)EX} (2) 
Symmetric structuring element is contained inside input set X as output of erosion is also translated 


as set of translation points. Figure 7 shows the example of dilation operation while Figure 8 represents 
erosion operation. 





(1) (il) 


Figure 7. Dilation operation (1) before operation (11) after operation 


Indonesian J Elec Eng & Comp Sci, Vol. 18, No. 1, April 2020 : 377 - 384 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 O 381 





‘@ (ii) 


Figure 8. Erosion operation (1) before operation (11) after operation 


Furthermore, the segmented disk image showing optic disc using Adaptive Thresholding method is 
secured. By using this method, optic disc can be carefully and clearly be extracted out. Optic disc is the most 
important characteristic in diagnosing glaucoma [17]. This owe to the fact that glaucoma is very hard to be 
roughly detected with only our eyes. It needs proper eye examination with suitable machine. Optic disc can 
show one of the obvious sign in diagnosing glaucoma. 

Lastly, the segmentation was made to identify optic cup as well. Optic cup is also very important 
characteristics in glaucoma detection. Its size can determine whether an eye is healthy or glaucomatous. 
From images obtained, it is proven that the optic cup of glaucomatous eye is bigger than optic cup of normal 
eye [18-20]. All the images for each process can be viewed in Table 1. 

Circular Hough Transform (CHT) method was commonly implemented as segmentation of optic 
disc [21]. By comparing to the proposed method, it can be seen that Adaptive Thresholding and Morphology 
method is more accurate as it also extract out optic cup and optic disc of retina fundus image. The CHT 
method only capable to show and illustrates countless circles as in shown in Figure 9 (1) while in Figure 9 
(11), one circle to detect optic disc was picked as maximum energy in the vicinity. 


Table 1. Images Obtained after Each Process 
RIM-ONE DRIONS-DB 
Normal eye Glaucomatous eye Normal eye Glaucomatous eye 


1 
_ Es 

e 

7 


a 
Ee 
Le 


the best approximate circle (maximum energy in the vicinity) [21] 





Original 
image 


Grayscale 
image 


Binary 
disc image 






Segmented 
disc image 


Binary cup 
image 


Segmented 
cup image 
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(1) (11) 


Figure 9. (1) Demonstrates the various circles (in red) detected by the Hough Transform (11) Shows 


2.3. Performance Metrics 

There are two parameters calculated for performance metrics for pre-processing. The Peak Signal to 
Noise Ratio (PSNR) and the Mean Square Error (MSE) are the two error metrics used to compare image 
compression quality. The calculation for PSNR and MSE is completed to compare the accuracy of the pre- 
processing method used [22-23]. 

Peak signal-to-noise ratio (PSNR) shows the ratio of power limit to the power of same image with 
noise. The PSNR calculate the peak signal-to-noise ratio between two images. It 1s normally regarded in 
logarithmic decibel scale (dB). However, it also can be calculated in percent (%). MAXI (maximum possible 
pixel value) of the image used 1s 255 because the image 1s uint8 [24]. 

As shown in (3) below is the formula for calculating PSNR and (4) is the formula to 
calculate MSE. 


PSNR =10 x{(log]_10 ((MAXI]‘2/MSE) (3) 


MSE = I/(mxn}>)_ G=0)4(m- 1)" G=0)(n- 1) [Image(i,j)-NoisyImage(1,})])|*2 ) (4) 


3. RESULTS AND ANALYSIS 

It can be concluded that for RIM-ONE database, the PSNR is 0.01891 and MSE is 65.62625. 
Meanwhile, for DRIONS-DB database, the PSNR 1s 64.0928 and the MSE is 0.02647. The results obtained 
clearly stated that the PSNR value obtained from after pre-processing is higher than the original image. The 
higher the PSNR value, the higher the level of audio signal than level of noise [25]. So, it 1s proven that the 
image quality after pre-processing is better than original image as noise was removed [26]. Table 2 indicates 
the PSNR value of original image before and after pre-processing was made. 


Table 2. The PSNR Value Before and after Pre-Processing 


RIM-ONE DRIONS-DB 
Image MSE PSNR Image MSE PSNR 
sample value (%) — value (%) sample from value (%) value (%) 
from database 
database 
N Im133 0.0158 66.1523 image_001 0.0366 62.4979 
O 
i Im152 0.0196 65.2033 image_002 0.0314 63.1555 
m Im160 0.0197 65.1949 image_003 0.0348 62.7135 
a Im161 0.0246 64.2238 image_031 0.0204 65.0375 
l Im169 0.0278 63.6834 image_062 0.0113 67.6144 
G Im013 0.0103 68.0137 image_006 0.0283 63.6138 
l 
a Im016 0.0120 67.3244 image_008 0.0275 63.7350 
u Im038 0.0151 66.3510 image_024 0.0239 64.3449 
c Im103 0.0309 63.2242 image_047 0.0246 64.2258 
O Im120 0.0133 66.8915 image_092 0.0259 63.9897 
m 
a 
Average (%) 0.01891 65.62625 0.02647 64.0928 


4. CONCLUSION 


As conclusion, the segmentation and detection of glaucoma is successfully completed. The methods 
used for preprocessing and segmentation works perfectly and are able to give satisfactory result. For future 
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work, the next process must be conducted. For instances, next process should involve calculating the size of 
optic cup and optic disc then proceed to classification into two common different types of glaucoma. Besides, 
a MATLAB GUI can also be developed in the future for a more efficient and accurate glaucoma detection 
system. 
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